Band structure of NiO revisited 
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The band structure of a strongly correlated semiconductor as NiO has been the object of much 
debate [PRL 103, 036404 (2009); PRL 102, 226401 (2009)]. Most authors, using computational 
techniques well beyond the simple density functional theory and the approximations GGA or LDA, 
claim that the band gap is about 4.0 eV and that the conduction band is of NiSd nature. Thus they 
seem to forget the results of electron energy-loss spectroscopy and inelastic x-ray scattering, both 
able to determine electronic transitions of only about 1.0 eV to an optically forbidden NiSd band. 
Further, a simple atomic calculation of the Ni"*""^ spin flip energy demonstrates that a NiSd band at 
4.0 eV is impossible. To set the issue straight, we calculated NiO with the very successful technique 
of PRB 78, 125116 (2008). It turns out that a band at 4.0 eV is optically accessible and made 
of excited atomic states, not Ni3d. Aside from that, we also found a narrow NiSd band at about 
1.0 eV. To confirm our procedures once again, we also calculated MnO and obtained the standard 
results of the good calculations as those cited above, and of experiment. 
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The band structure of NiO is being the object of many calculations by different techniques. To quote just a few, 
Physical Review Letters recently published the works of Engel and Schmid [l| on the Exact-Exchange technique, 
and of Tran and Blaha on a novel semilocal exchange-correlation potential. Physical Review published the work 
of Rodl et al 3 on a GGA calculation. The Journal of Physics: Condensed Matter published a self-interaction 
(SIC) calculation [|. The many authors claim that NiO presents an important challenge for density functional theory 
because it is a strongly correlated system. The most simple LSDA and GGA calculations are reported to be failures. 
Thus, along the years NiO is being calculated with many sophisticate techniques. Aside from the works already 
quoted, we can cite some examples: Deng et al used a variational LDA formulation 0, Fuchs et al0| with their 
generalized Kohn-Sham method, Kobayashi et al [1] that used a GW started from LSDA+U, Eder Q^that used a 
Variational Cluster Approximation, Miura et al Kuncs ct al 11 1, and Ren et al []3 using different LDA+DFMT, 



Tran et al [l^ with a hybrid exchange-correlation, mixing Hartree-Fock, Han et al |l4l| with LDA+U, Li ct al [If 
and Kotani and van Schilfgaarde [ll] using a self-consistent GW, Korotin et al [l3| with a Wannier function method. 

Of course it is impossible to review the many papers. We want only to call attention that many of these calculations 
are wrong when they set the band gap at about 4.3 eV, following the experimental BIS result of Sawatzky and Allen 
[isj . claiming that the conduction band at that energy is a Ni-d band. The location of the conduction Ni-d bands has 
been determined by Miiller and Hiifner [Toj using electron energy-loss spectroscopy, and by Huotari et al po| using 
inelastic x-ray scattering. The band at 4.3 eV above valence is certainly not Ni-d because these start at a much lower 
energy, about 1.0 eV. It is simple to show the approximate position of the Ni-d conduction bands by calculating the 
following atomic energy differences 



or the Ni++ ionic difference 



7V^[4s2(3d t)^(3d i)^] - NtiAs^Sd t)'(3d if 



m[{3d t)^(3d If] - Ni[{3d T)'(3d If 



These differences are the energies required to flip a spin, which correspond to the excitation of an electron to the 
antiferromagnctic crystal conduction band. Using the PBE exchange-correlation the difference is 0.99 eV for 
the atom, and 1.03 eV for the ion. Using the LDA exchange-correlation of Cepperley- Alder Perdew-Zunger [2l[ the 
differences are only 0.05 eV smaller. In any case, the spin-flip en ergy is similar to the threshold to the first energy 
band, as found experimentally [l^, and the 4.3 eV BIS band [lj| could never be a Ni-d band. Compared with 



'Electronic address; |guimaOO@gmail.com| 

t Electronic address: lkteles@it a.brl 

t Electronic address: mmarques@ita.br| 



2 



MnO, also an antiferromagnctic oxide. NiO is very different. In tlie case of MnO, the spin-flip energy 



is 3.55 eV, similar to the band gap, so that a Mn-d conduction band is possible. 

Next we present the results of our calculation for NiO. Our calculation was all-electron with wavefunctions expanded 
in LAPW-l-lo [13] ■ We avoided the pseudopotential calculations for two reasons. First, in the case of NiO, the Ni-d 
functions are not easily expanded with plane waves. Second, for magnetic systems, the exchange-correlation potential 
may depend much on the way the core density is treated [23|, because linearization is unacceptable. So it is not just a 
problem of finding the good pseudopotential, but the good core density must also be found. NiO is an antifcrromagnet 
with the rock salt structure. The Ni sublattice is antiferromagnetically ordered as an Lli (CuPt) alloy. It has a very 
small trigonal distortion along the body diagonal. Along this direction, the planes alternate spins, up and down. We 
used the experimental lattice parameters and did not optimize geometry, which is the best procedure to calculate the 
elementary excitation spectrum. 

In Fig. [Ija) we show the valence and conduction bands of the Kohn-Sham eigenvalues calculated with the PBE 
exchange-correlation Contrary to the opinion of many authors, we do not see anything wrong with those bands. 
The Kohn-Sham eigenvalues should not compare well with the excitation spectrum, the band gaps being notoriously 
smaller. Unfortunately there seems to be no procedure to simultaneously calculate the total energy, which is the object 
of Density Functional Theory, and the excitation spectrum. The spectrum calculation should be made following a 
very good GW [l^l that uses the self-energy operator instead of the exchange-correlation potential. Though not 
correct, the Kohn-Sham energy spectrum usually gives an indication of how the correct spectrum will be. As it stands 
in Fig. [n^a) NiO is an insulator, the first conduction bands arc about 1 eV above the valence band and is made of 
minority Ni d wavefunctions. The identification of the many valence states as Ni d majority and minority spins follows 
the simple argument: atomic Ni has 2 electrons 4s and 8 electrons 3d. In NiO, the nickel looses the two electrons 4s to 
the Oxygen, there completing its 2p shell. Due to the cubic crystal field, the five 3d states of Nickel splits into ^2 with 
degeneracy 3 and e with degeneracy 2. Due to the trigonal distortion, t2 is further split into states with degeneracies 
2 and 1. Of course these identifications can only be seen at the zone center F. Thus among the 8 remaining Ni valence 
electrons, 3 occupy t2 with majority spin, 2 occupy e with majority spin, and 3 occupy t2 with a minority spin. The 
state e with minority spin is empty and defines the first conduction band. 

The first conduction bands are really acceptor bands. An electron added to them would change the charge state of 
the Ni ion, according to 



These acceptor bands conduct very poorly and only by hop ping from one Ni ion to another with the same spin. 
They are responsible for the dd excitation spectra of [l9[ and [20| . at energies very different from the Sawatzky and 
Allen gap of 4.3 eV. The band gap to these acceptor bands is only 1.00 eV, as calculated by GGA [1], or 0.51 eV as 
calculated with LDA [2^. We must mention that one sees this gap even smaller in the literature due to the use of 
plane wave bases. Anyway, though Fig. [Ija) gives the correct description of NiO, it misses the Sawatzky and Allen 
gap altogether. The bands of NiO, described by Fig.[TJa), and the results of Sawatzky and Allen remain incompatible. 

Instead of a very precise GW calculation, Ferreira et al [1^ developed a simple and very successful procedure, with 
no adjustable parameters, to calculate the excitation energy spectrum. The procedure is inspired in the old transition 
state technique for atoms, shown to be equivalent to the inclusion of the self-energy of the quasi-particle. Here it 
must be recalled that the "self-energy" being calculated is not that of Hedin [23j, but the classical self-energy of 
Electrostatics 



to which a minor exchange-correlation term is added. Thus the self-energy operator of GW is expected to be equivalent 
to the exchange correlation potential plus the self-energy potential defined in p6j , and plus an imaginary term that 
has eluded all calculations. As with the best GW calculations, the Ferreira et al method, named LDA-1/2, produces 
very good band gaps and effective masses, but at a very small computational price. Lately we have been also using a 
GGA-1/2 method, where the exchange-correlation is GGA, instead of LDA, and the -1/2 is to remind that we remove 
1/2 electron, as in the transition state technique. The self-energy of [2^ has a very simple physical meaning: it is the 
work necessary to join all the charge of the quasi-particle, dispersed in a Bloch state, into a localized wavefunction. 
As such, this energy has to be subtracted from the Kohn-Sham eigenvalue. The self-energy potential is calculated in 
the atom, and then trimmed so that it does not extend to neighboring atoms. The trimming is made by means of a 
cutting function with a parameter "CUT" which is determined variationally by making the band gap extreme. In the 
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case of NiO we used CUT{NiM) = 1.97a.u. and CUT{02p) = 2Ala.u. to make the "4.3 band gap" maximum. The 
maxima are not sharp and neighboring values of CUT could be used as well. 

The energy spectrum calculated with the inclusion of self-energies is represented in Fig.[ljb). The results obtained 
with GGA does differ from those from LDA. Aside from the relative upwards motion of the "4.3 band" , there is not 
much difference from the results of a standard all-electron calculation (Fig. [TJa)). It is clear why the "acceptor band" 
at ~ leV does not move with respect to the valence bands: Both are Ni3d bands and have the same self-energy. On 
the other hand, these acceptor bands do not exist in the case of MnO. In this latter case, using our technique, we 
obtained a gap of 4.18 eV (GGA) and 3.98 eV (LDA), matching the results of other calculations and experiment, 
and no surprises were obtained [23| . 

To conclude we want to say that standard LDA or GGA calculations seem reliable even for a strongly correlated 
material as NiO. But it is important to remember that DF will never give the excitation spectrum. An attempt 
to force the obtaining of such spectrum may generate errors and misinterpretation of the calculated results, as the 
impossible 4.3 eV Ni d band. 
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FIG. 1: (a) NiO energy bands calculated with GGA. (b) NiO energy bands calculated with GGA-1/2, and in parenthesis, 
calculated with LDA-1/2. 
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